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^ ■ Abstract 

The linear /c^-factorization is part and parcel of the perturbative QCD description of high 
Q^' energy hard processes off free nucleons. In the case of heavy nuclear targets the very concept of 



nuclear parton density becomes ill-defined as exemplified by the recent derivation [2] of nonlin- 
ear nuclear ^-factorization for forward dijet production in deep inelastic scattering off nuclei. 



Here we report a derivation of the related breaking of /^-factorization for single-jet processes. 
We present a general formalism and apply it to several cases of practical interest: open charm 
and quark & gluon jet production in the central to beam fragmentation region of ^y*p,^*A,pp 
and pA collisions. We show how the pattern of fcj_-factorization breaking and the nature and 
number of exchanged nuclear pomerons do change within the phase space of produced quark 
and gluon jets. As an application of the nonlinear fc^-factorization we discuss the Cronin effect. 
Our results are also applicable to the p^-dependence of the Landau-Pomeranchuk-Migdal effect 
for, and nuclear quenching of, jets produced in the proton hemisphere of pA collisions. 

PACS: 13.87.-a, 13.85.-t, 12.38.Bx, 11.80.La 
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I. INTRODUCTION 



The familiar perturbative QCD (pQCD) factorization theorems tell that the hard 
scattering observables are linear functional (convolutions) of the appropriate parton 
densities in the projectile and target [1] and are based on an implicit assumption that in 
the low-parton density limit the unitarity constraint can be ignored. They play, though, 
a fundamental role in interaction with strongly absorbing heavy nuclei which puts into 
question the very concept of a nuclear parton density. Indeed, a consistent analysis 
of forward hard dijet production in deep inelastic scattering (DIS) off nuclei revealed 
a striking breaking of /c^-factorization [2]: Namely, starting with the nuclear-shadowed 
nuclear DIS cross section, one can define the collective nuclear unintegrated gluon density 
such that the familiar linear /c^-factorization (see e.g. the recent reviews [3]) would hold 
for the forward single-quark spectrum but fail for the quark-antiquark two-jet spectrum 
which turns out to be a highly nonlinear functional of the collective nuclear gluon density. 
This finding casts a shadow on numerous extensions to nuclear collisions of the parton 
model formalism developed for nucleon-nucleon collisions and calls upon for revisiting 
other hard processes. The ultimate goal is the consistent description of nuclear effects 
from DIS to ultrarelativistic hadron-nucleus to (the initial stage of) nucleus-nucleus 
collisions, for instance at RHIC and LHC. 

The purpose of this communication is a demonstration that the pattern of k±- 
factorization breaking discovered in [2] holds also for a broad class of single-jet spectra. 
Here we recall that nuclear shadowing in DIS sets in at the Bjorken variable 

1 

x<x A = - , (1) 

R A m N 

when the coherency over the thickness of the nucleus holds for the qq Fock states of the 
virtual photon ([4, 5], for the color dipole phenomenology of the experimental data on 
nuclear shadowing see [6]). Here Ra is the radius of the target nucleus of mass number 
A and m N is the proton mass. The analysis [2] focused on the excitation of quark and 
antiquark jets, 7* — > qq, in the photon fragmentation region at x < Xa- The subject of 



2 



this communication is excitation of open charm and jets, 

g* -> QQ, g* -> 00, q* -»• (2) 

in hadron-nucleon/nucleus scattering and/or DIS off nucleons and nuclei at rr <C xa- At 
RHIC that means the large (pseudo)rapidity region of proton-proton and proton-nucleus 
collisions. Our starting point is the color dipole multiple scattering theory formulation 
of the open charm production in pA collisions [7] which we extend from the total to 
differential cross section and to arbitrary color state of the projectile parton. We treat 
the jet production at the pQCD parton level, our approach to the jet spectra is in many 
respects similar to Zakharov's lightcone description of the QCD Landau-Pomeranchuk- 
Migdal (LPM) effect for the limit of thin targets [8, 9]. The results obtained below for the 
open charm and quark and gluon jets make manifest a change from linear to nonlinear 
A;_i_-factorization as the target changes from a free nucleon to heavy nucleus. 

In the production of dijets in DIS [2] and irA collisions [10] the fragmenting particle 
was a color-singlet one. Here we consider the gN, gA, qN, qA collisions in which the beam 
partons are colored. In general the scattering of colored particles is plagued by infrared 
problems, but the excitation cross sections can be calculated as combinations of infrared- 
safe total cross sections for interaction with the target of certain color singlet multiparton 
states [2, 7, 8]. We discuss to some detail the so-called Cronin effect [11] - the nuclear 
antishadowing at moderate to large transverse momentum p of jets. Decomposing in 
the spirit of the Leading Logp 2 (LLp 2 ) approximation the single-jet cross section into 
the so-called direct and resolved interactions of the incident parton we demonstrate 
that the Cronin effect is a salient feature of the resolved interactions and identify its 
origin to be the antishadowing nuclear higher-twist component of the collective gluon 
density of a nucleus derived in [12, 13]. Regarding the potential applications to pA 
collisions at RHIC we note that the /cj_-factorization phenomenology of DIS [14] suggests 
fairly moderate Q A ~ 1 GeV at x ~ 10~ 2 for the heaviest nuclei [2], but we find 
that interesting nuclear effects extend to the pQCD domain of jet transverse momenta 
way beyond Qa- Qualitatively similar results would be expected in a regime of strong 

3 



absorptive corrections also for hadronic targets other than nuclei, the nucleus represents 
a theoretical laboratory, in which saturation effects are endowed with a large parameter, 
the nuclear opacity, which allows for their systematic account. 

One more point needs an emphasis. The partonic cross sections evaluated within 
the often used collinear approximation are subject to large smearing corrections for the 
so-called intrinsic transverse momentum (p T ) of partons, see [15]. In contrast to that we 
treat the pQCD partonic subprocesses with full allowance for the transverse momenta 
of all interacting partons and no extraneous (pr)-smearing is needed. The discussion of 
nonperturbative effects such as nuclear modifications of the hadronization of jets ([16] 
and references therein) or of the recombination of partons into high-p^ hadrons ([17] and 
references therein) goes beyond the scope of the present study. 

All our results are for the distribution of jets in both the transverse and longitudinal 
momentum, the latter is parameterized in terms of the fraction, z, of the longitudinal 
momentum of the incident parton a carried by the observed parton b. As such, they solve 
the LPM problem for finite transverse momentum p, while the early works on the LPM 
effect focused on the p-integrated longitudinal momentum spectrum. The p-dependent 
LPM effect is important for a quantitative pQCD treatment of the nuclear quenching 
of forward high-p hadrons observed experimentally by the BRAHMS collaboration [18]. 
The numerical results for the LPM effect for forward jets will be reported elsewhere, 
here we focus on the derivation and implications of the nonlinear /c_i_-factorization for 
the transverse momentum distributions. 

The presentation of the main material is organized as follows. Our starting point is the 
color-dipole S'-matrix representation for the two-body spectrum from the fragmentation 
a — > be. We demonstrate how the single-parton cross spectrum can be calculated in 
terms of the S'-matrices for special 2- and 3-body color singlet parton states (some 
technical details are given in Appendix A). The color-dipole master formulas for the 
single-parton spectra in interactions with free-nucleon and heavy-nucleus targets are 
reported in section III. They form the basis of the /c^-factorization representation for 
the spectra of different partons and in section IV we start with the linear, and in section V 
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with nonlinear, /c^-factorization for open charm production off free nucleons and nuclei, 
respectively. The subject of section VI is the Cronin effect in open charm production 
off nuclei. We demonstrate that in the general case the nuclear spectrum is a quadratic 
functional of the collective nuclear gluon density and illustrate this property in terms 
of the Kancheli- Mueller diagrams for inclusive cross sections [19]. The Cronin effect is 
shown to be a generic feature of the so-called resolved gluon interactions and its origin is 
attributed to the nuclear antishadowing property of the collective nuclear gluon density 
derived in [12]. We show that despite the relatively small saturation scale Qa the Cronin 
effect extends to perturbatively large transverse momenta. In section VII we apply our 
master formula to the transverse momentum spectrum of gluons from the fragmentation 
9* - ► 99- 111 this case the nuclear cross section is a cubic functional of the collective 
nuclear gluon density. The subject of section VIII is the production of quark and gluon 
jets in the fragmentation q* — > qg. The relevance of our results to the p-dependent LPM 
effect is pointed out in sections VII and VIII. The main results are summarized and 
discussed in section IX. Our principal conclusion is that the /c^-factorization breaking is 
a universal feature of the pQCD description of hard processes in a nuclear environment. 
We comment on earlier works on single quark and gluon jet production [20-34], some 
of which invoke the Color Glass Condensate approach ([35] and references therein). A 
convenient technique of derivation of the n-parton interaction S-matrices is presented in 
Appendix A. The rules for the derivation of the wavefunctions of two-parton states in 
terms of the familiar splitting functions are given in Appendix B. 

II. INCLUSIVE PRODUCTION AS EXCITATION OF FOCK STATES OF THE 
BEAM PARTON 

Previously, the breaking of linear /^-factorization had been demonstrated for the ex- 
citation of dijets, 7* — > qq, in DIS at x < xa [2]. Specifically, the two-parton transverse 
momentum spectrum was shown to be a highly nonlinear functional of the collective 
nuclear gluon density. On the other hand, upon integration over the transverse mo- 
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menta of the antiquark jet, the single-quark spectrum was found to fulfill the linear 
/c^-factorization. In this paper we extend the analysis of [2] to the case of colored beam 
partons, and confirm the breaking of A;_i_-factorization as a generic feature of the pQCD 
description of hard processes in a nuclear environment. 

To the lowest order in pQCD the underlying subprocess for the mid- to large-rapidity 
open charm production in proton-proton collisions is the fusion of gluons from the beam 
and target, gg — > cc; for the gluon jet production one would consider a collision of a 
beam quark with a gluon from the target, qg — > qg; at higher energies the relevant 
subprocess for gluon production will be gg — > gg. All the above processes are of the 
general form ag — > be and, from the laboratory frame standpoint, can be viewed as an 
excitation of the perturbative \bc) Fock state of the physical projectile \a) by one-gluon 
exchange with the target nucleon. In the case of a nuclear target one has to deal with 
multiple gluon exchanges which are enhanced by a large thickness of the target nucleus. 
A general treatment of multiple gluon exchanges in nuclear targets has been developed 
in [2, 7], its extension to the a — > be transitions can be described as follows: 

Partons with energy E a and transverse momentum p a such that 

2m N E a 1 
Vl X A 

propagate along straight-line, fixed-impact-parameter, trajectories and interact coher- 
ently with the nucleus, which is behind the powerful color dipole formalism [5, 36, 37]. 
The target frame rapidity structure of the considered a — > be excitation is shown in 
fig. 1. The beam parton has a rapidity r) a > t\a — logl/au, the final state partons too 
have rapidities r] b c > tja- In this paper we focus on the lowest order excitation processes 
a — > be without production of more secondary partons in the rapidity span between i] b 
and r) c . 

To the lowest order in the perturbative transition a — > be the Fock state expansion 
for the physical state \a) p h ys reads 

\a) P hys = \a)o + ^(z b ,r)\bc) (4) 
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FIG. 1: (Color online) The rapidity structure of the excitation a — > be: (a) excitation of mid- 
to large-rapidity open charm g — > QQ, (b) radiation of gluons by quarks q — ► qg, (c) radiation 
of gluons by gluons g ^ gg- 



where ^{z b) r) is the probability amplitude to find the be system with the separation r 
in the two-dimensional impact parameter space, the subscript "0" refers to bare partons. 
The perturbative coupling of the a — > be transition is reabsorbed into the lightcone wave 
function ^>(zi,r), and we also omitted a wave function renormalization factor, which is 
of no relevance for the inelastic excitation to the perturbative order discussed here. 
If b a = b is the impact parameter of the projectile a, then 

b b = b + z c r, b c = 6 - z b r , (5) 

see Fig. 2. Here z byC stand for the fraction the lightcone momentum of the projectile 
a carried by the partons b and c. The virtuality of the incident parton a is given by 
Qa = tfi-, where k a is the transverse momentum of the parton a in the incident proton 
(Fig. 1). For the sake of simplicity we take the collision axis along the momentum of 
parton a. The transformation between the transverse momenta in the a-target and p- 
target reference frames is trivial, pfl = p^l + z c ^k a: below we cite all the spectra in the 
a-target collision frame. We shall speak of the produced parton - or jet originating from 
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FIG. 2: (Color online) The color dipole structure (lhs) of the generic 3-parton state and (rhs) of 
the 3-parton state entering the color-dipole description of fragmentation a — > be with formation 
of the be dipole of size r. 



this parton - as the leading one if it carries a fraction of the beam lightcone momentum 
z — *■ 1 and a slow one if it is produced with 2 C 1. We speak of the produced parton 
as a hard one if it is produced with large transverse momentum and a soft one if its 
transverse momentum is small compared to the so-called nuclear saturation scale to be 
denned below. 

By the conservation of impact parameters, the action of the S-matrix on \a) p f iys takes 
a simple form 

S\a) phys = S a (b)\a) + S b (b b )S c (b c )^(z b ,r)\bc} 

= S a (b)\a) phys + [S b (b b )S c (b c ) - S a (b)]*(z b ,r)\bc) . (6) 

In the last line we explicitly decomposed the final state into the elastically scattered 
\o)phys and the excited state \bc). The two terms in the latter describe a scattering 
on the target of the be system formed way in front of the target and the transition 
a —>■ be after the interaction of the state \a)o with the target, as illustrated in fig. 3. The 
contribution from transitions a — > be inside the target nucleus vanishes in the high-energy 
limit 

x ~ wh b - Xa (7) 

where M bc is the invariant mass of the excited be system, Q\ is the virtuality of the 
incident parton a, W is the total energy in the beam-target nucleon collision center of 
mass frame and xb is the fraction of the beam energy carried by the beam-parton a. 
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a) b) c) 

FIG. 3: (Color online) Typical contribution to the excitation amplitude for gA — ► qqX, with 
multiple color excitations of the nucleus. The amplitude receives contributions from processes 
that involve interactions with the nucleus after and before the virtual decay which interfere 
destructively. 

The probability amplitude for the two-jet spectrum is given by the Fourier transform 

J d 2 b b d 2 b c eM-KP b b b + P c b c )][S b (b b )S c (b c ) - S a {b)}*{z b ,r) (8) 

The differential cross section is proportional to the modulus squared of (8) and one 
encounters the contributions containing 

s£\b',b) = st(b')s a (b) 

S®(b',b h b c ) = Sl(b')S b (b b )S c (b c ), 

sQ a (b,b' b ,K) = st(b' b )sl(K)s a (b) 

S^ b K b' c , b b , b c ) = St(b' b )Sl(b' c )S c (b c )S b (b b ) . (9) 

Here we suppressed the matrix elements of over the target nucleon, for details see 
[2] . In the calculation of the inclusive cross sections one averages over the color states of 
the beam parton a, sums over color states of final state partons b, c and takes the matrix 
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FIG. 4: (Color online) The S-matrix structure of the two-body density matrix for excitation 
a — ► 6c. 



products of and S with respect to the relevant color indices entering S^ n \ Some of 
the technicalities of the derivation of are presented in the Appendix A, here we cite 
the resulting two-jet cross section: 

dzd 2 p b d 2 p c (27r) 4 J 

x exp[ip b (b b - b' b ) + i Pc (b c - b' c )]V{z, b b - b c )V*(z, b' b - b' c ) 

{S^K b' c , b>, b c ) + S^(b\ b) - S&(b, b' b , b' c ) - S%l(b', b b , b c )} (10) 

The crucial point is that the hermitian conjugate can be viewed as the S'-matrix for 
an antiparton [2, 7, 38]. As a result, S^(b',b) is an S'-matrix for interaction with the 
target of the ad state in which the antiparton a propagates at the impact parameter b'. 
The averaging over the color states of the beam parton a amounts to taking the color 
singlet ad state. The S^ bc and Sr_ , describe interaction of the color-singlet abc and bccb 
states, respectively. This is shown schematically in Fig. 4. 

The and are readily calculated in terms of the 2-parton and 3-parton dipole 
cross sections [5, 7, 37], general rules for multiple scattering theory calculation of the 
coupled-channel are found in [2] and need not be repeated here. 

The results for the two-jet cross sections will be presented elsewhere, here we focus on 
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the single-jet problem. Integration over the transverse momentum p c of the jet c gives 

b c = K, (11) 

b-b' = z b (r-r'), b h -b' b = r-r', b'-b c = z h r'. (12) 

In the generic case the 3-parton state abc has the dipole structure shown in Fig. 2 and 
the dipole cross section a 3 abc(p, r )- In the considered problem of the a — > be excitation 
p = z b r'. The case of S^ cb deserves special scrutiny. In the general case, is a 
multichannel operator, an example is found in [2] . Because in the single-particle spectrum 
b' c = b c , the unitarity relation 

Sl(b c )S c (b c ) = 1 (13) 

leads to a fundamental simplification of S^, b (b b ,b' c ,b b ,b c ). Specifically, in a somewhat 
symbolic form 

S£> cb (b' b , b' c , b b , b c ) = Sl(b' b )Sl(b c )S c (b c )S b (b b ) = Sl(b' b )S b (b b ) = S£\b' b , b b ) , (14) 

i.e., the effect of interactions of the spectator parton c vanishes upon the summation 
over all its color states and integration over all its transverse momenta [7]. The coupled- 
channel operator then takes a single-channel form, some technical points behind 
this derivation are clarified in Appendix A. The only trace of the observed parton b 
having been produced in the fragmentation a — > bc is in the density matrix ^(z b , b b — 
b c )^*(z b ,b' b — b' c ) which defines the transverse momentum distribution of the parton b 
in the beam parton a and the partition z b , z c = (1 — z b ), of the longitudinal momentum 
between the final state partons. It will be more convenient to regard S^ b \b b ,b b ) as a 
function of the transverse size of the bb dipole b' b — b b , and the impact parameter B of the 
dipole-nucleus interaction, and we shall from now on write Sj?\b' b , b b ) — > 5^(5, b' b — b b ). 

The observation (14) can be generalized further to the cancellations of beam spectator 
interactions. The detailed discussion will be reported elsewhere, here we just quote the 
simple example of open heavy flavour excitation. The above expounded formalism for 
the fragmentation g — > QQ can readily be generalized to excitation of heavy flavour from 
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the incident quark q* — > g'QQ which proceeds via q* — > g'g, 5 — > QQ- The simplest 
case of direct importance for mid-rapidity heavy flavour production in pp,pA collisions is 
when the produced QQ pair carries a small fraction of the incident quark's momentum, 
i.e., the underlying partonic subprocess is a collision of a slow gluon from the beam with 
a slow gluon from the target, gg — > QQ. The generalization of eqs. (11), (12) to this 
case shows that the beam quark q* and the spectator quark q' would propagate at the 
same impact parameter b q . One readily finds that the amplitude of the q* — > q'QQ 
will be given by eq. (6) times the extra factor of S q (b q ). Upon the integration over the 
transverse momentum of the spectator quark one gets the equality of the corresponding 
impact parameters in the multiparton ^-matrix and its hermitian conjugate, eq. (11), 
and upon the application of the unitarity relation (13) the effect of spectator interac- 
tions cancels out in complete analogy to eq. (14). The resulting Sj?\B,b' b — b b ) only 
depends on the dipole parameter b' b — b b and the overall impact parameter, B, of the 
multiparton system. What will be left of the beam quark q* is the wave function of its 
q'g Fock state which defines the longitudinal and transverse momentum density in the 
beam of the gluon g* which excites into the QQ state. The evident assumption behind 
the spectator interaction cancellations is that the coherency condition is fulfilled for all 
spectator partons; because the spectators are faster than the observed parton, this weak 
condition does not impose any stringent constraints beyond (7). 



III. COLOR DIPOLE MASTER FORMULA FOR FREE NUCLEONS AND 
NUCLEI 



For the free nucleon target integration over the impact parameter B gives the corre- 
sponding total cross sections, for instance, 



2 J SB [l - SH\B, b' b - b b )} = a 2 ~ bb {r - r') 



2 d 2 B 



1 - S®(B, b' - b c , b b - b c ) = a 3 - abc (z b r', r) . 



(15) 
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Then the single-jet spectrum for the free nucleon target can be cast in the form 
da N (a* -> b(p b )c) 1 f „ , 2 , 



J d 2 rd 2 r' exp[ip b (r - r')]V(z b ,r)V*(z b ,r') 



dz b d 2 p b 2(2tt) 2 
x{a 3! abc(z b r\r) + a 3 ^(z b r,r') - a 2;bb (r - r') - a 2 ^(z b (r - r'))} (16) 

In the generic dbc system there could be more than one scheme of coupling to a color 

(3) 

singlet state, in which case S^ bc would be a coupled-channel operator, but that is not the 
case with gluons and quarks. The coupled-channel operator S^ ch also takes the single- 
channel form (14). As usual, we apply multiple scattering theory treating the nucleus as 
a dilute gas of color-singlet nucleons, and upon the summation over all the nuclear final 
states and application of the closure relation the nuclear matrix elements of take 
the familiar Glauber-Gribov form [39, 40]. The resulting nuclear single-jet spectrum per 
unit area in the impact-parameter plane equals 

x [T A [b, (T 3 ,abc(z b r' , r)} + T A [b, a 3 ^(z b r, r')} 

-r A [b,a 2tbb (r - r')\ - T A [b,a 2 , a - a (z b (r - r'))]} . (17) 
We can use the Glauber-Gribov multiple scattering theory result 

T A [b,a] = l-exp[-^T(b)], (18) 

/oo 
dr z n A (b,r z ) is the optical thickness of the nucleus, and the nuclear 
-oo 

density n A (b,r z ) is normalized according to J d 3 fn A (b,r z ) = J d 2 bT(b) = A. 

Within the color-dipole formalism the x-dependence of the DIS structure function is 
driven by the x-dependence of the color-dipole cross section <J q q(x, r) which is governed by 
the color-dipole form of the Balitskii-Fadin-Kuraev-Lipatov (BFKL) equation [37, 41, 
42]. For x > x A DIS off a nucleus amounts to the sum of incoherent scatterings off bound 
nucleons. The onset of nuclear coherence effects at x < x A must be treated with the 
color-dipole cross section a qq (x A) r) which has been BFKL-evolved down to x — x A . We 
recall that for heavy nuclei well evolved shadowing sets in at x A w 1CT 2 . The considered 
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lowest order excitation processes a — > be with rapidities rj^ c > t\a such that no more 
secondary partons are produced in the rapidity span between rjt, and r\ c or between i] b c 
and t]a must also be treated in terms of a qq (xA, r). The values of x^ c attainable at 
RHIC are such that jet production in the proton hemisphere of pA collisions at RHIC 
is dominated by precisely these lowest-order processes. The allowance for higher-order 
processes would amount to the x-evolution of the /c_|_-factorization breaking and will be 
addressed to elsewhere. 



IV. /c ± -FACTORIZATION FOR OPEN CHARM PRODUCTION OFF FREE 
NUCLEONS 

Here the underlying pQCD excitation is g ^ QQ, i.e., a = g,b = Q,c = Q. The 
process is symmetric under Q <-> Q and for the sake of brevity we we put z = zq and 
p = pg. The relevant 2-parton and 3-parton dipole cross sections equal [7, 37] 

a 3 (x, zr', r) = [a qq (x, zr' - r) + a qq -(x, zr') - a qq -(x, r)\ + a qq (x, r) , 
C 

a gg (x,z(r-r')) = -± a qq (x, z(r - r')) , (19) 
N 2 - 1 

where Ca = N c , Cf = ^ AT — are the familiar quadratic Casimirs, and [43] 



Vqqfa r ) = J d 2 Kf(x, k)[1 - exp(inr)] , 



/(I , K) = !^M i . Hx y } . (20, 

iv c K 

<r qq (x, r) is the dipole cross section for the qq color dipole and 

^) = ^4> (21) 
a log k 

is the unintegrated gluon density in the target nucleon (our normalization of f(x A , k) and 
of the related Weizsacker- Williams gluon density 0(6, xa, k) = (f) WW (b, xa, k), see below, 
are slightly different from the ones used in [2]). In the momentum space calculations the 
running coupling as(r) must be taken at an appropriate hard scale. 
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Applying the technique developed in [2], upon the relevant Fourier transformations 
one readily finds 

2(2n) 2 da N (g* -> QQ) 



dzd 2 p 



= J d 2 Kf(x A ,K) 



N 2 

i m-l (l^^' p ) -V(z,p + zk)\ 2 + \V(z,p + n) -V(z, P + zk) 

iv c 

+ J^— [ \y(z,P)-y(z,P + K)\ 2 } (22) 

= J d 2 Kf{x A ,K)^^{\^{z,p)-^(z,p + ZK)\ 2 +\^{z,p + K)-^{z,p + ZK,)\ 2 

-\V(z,p) - V(z,p + k)\ 2 } + \V(z,p) - V(z,p + k)\ 2 ^J . (23) 

Here we show explicitly that the target gluon density f(x A , k) enters with the boundary 
condition argument x A . The expression \^(z,p) — ^(z : p + zk)\ 2 for the QQ Fock state 
of the gluon is obtained from its counterpart for the transverse photon cited in [2] by 
the substitutions Q 2 — > Q\ and N c a em (?Q — > Tpas, see Appendix B. 

Here we emphasize that Eqns. (22,23) account fully for the transverse momenta of all 
involved partons, no further smearing of the obtained spectra for the intrinsic transverse 
momentum of partons is needed, furthermore, it would be illegitimate. Early discussions 
of charm production in hadron-nucleon collisions in the framework of /c^-factorization 
can be found in [44], for reviews and a guide to recent literature we refer to [3]. The rather 
straightforward representation (22,23) of the heavy quark cross section in terms of light- 
cone wavefunctions for the transition g* — * QQ is new, as are the related expressions for 
quark/gluon production to be found below. 

Some comments on this simple formula are in order. First, although we spoke of the 
charm production, as soon as the transverse momentum p is perturbatively large eq. (23) 
will fully be applicable to light quark jets and there is an obvious connection between 
(23) and the NLL real correction to the BFKL kernel from light qq pair production, 
but dwelling into that is beyond the scope of the present communication. Second, we 
presented our result for two different groupings of terms. While in the first grouping 
(22) the result is given as a manifestly positive sum of squares, in the second one (23) we 
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decomposed the cross section into 'nonabelian' , oc Ca/2Cf, and 'abelian' pieces. Putting 
Ca — > one would switch off the nonabelian coupling of gluons to gluons and obtain 
indeed the same formula as for the single-quark spectrum from the 7* — > qq transition, 
i.e., for inclusive DIS [2]. Third, note that the same 'abelian' limit is obtained at the 
boundaries z = and 1 — 2 = 0. In the both cases that can be traced to the fact either 
the spectator quark for z = 0, or the observed antiquark for z = 1 would propagate 
at precisely the same impact parameter as the parent gluon so that the non-abelian 
color structure of the gQQ is simply not resolved. Still the physics of the two cases is 
different: in the limit of z = the 'initial state' gg dipole interaction effectively drops 
out, whereas at z = 1 there are subtle cancellations between the gg dipole cross section 
and the octet-octet dipole components of the 3-parton cross sections. As we shall see 
below, the latter cancellations will be upset for nuclear targets. Fourth, in the large- N c 
limit Ca = 2Cp and for later reference we quote 

2(2ir) 2 da N (g* -> QQ) 

dzd 2 p Wc>i 

= j d 2 Kf(x A ,K) {\V(z,p) -V(z,p + zk)\ 2 + \^(z,p + K) -V(z,p + zk)\ 2 } .(24) 

Based on the discussion of spectator effects in section III and anticipating the dis- 
cussion in section VI, we note that the open-charm production can be presented in a 
manifestly beam-target symmetric form. One can multiply the above derived cross sec- 
tion d<jq ~* QQ \z,xa,p) by the unintegrated flux of gluons of transverse momentum k a 
and lightcone momentum fraction xb in the beam proton, T{xb, k a ), take the wave func- 
tions of Appendix B at Q 2 = k 2 , and calculate the charm spectrum for pp,pA collisions 
as 

1 f * 
da N)A {z,x A ,p) = — / d 2 k a J r (x B ,k a )d(T%A QQ \ z > x A,P + zk a ) . (25) 
Xb J 

The usual kinematical constraint, xbXaW 2 = —. ^ , is understood. The fully inclu- 

z(l-z) 

sive single-parton cross section is obtained from (25) upon integration over xb keeping 
fixed the longitudinal momentum of the observed parton zxb- We stick to z because (i) 
the transverse momentum spectra exhibit a highly nontrivial dependence on z and (ii) 
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c) d) 



FIG. 5: (Color online) (a) The Kancheli-Mueller diagram for the two-parton inclusive cross 
section in proton-proton collisions, (b) its representation in terms of the pomeron exchange 
and (c) transition to the optical theorem for the single-parton inclusive cross section (d) upon 
the integration over the phase space of the parton c. 

understanding the z-dependence is crucial for the LPM effect. 

Within the fcj_-factorization the microscopic QCD description of the pomeron ex- 
change is furnished by the unintegrated gluon densities. To this end one can look at 
the open-charm cross section (25) through the prism of the Kancheli-Mueller diagrams 
[19] for inclusive spectra as depicted in a somewhat symbolic way in figure (5). The 
cbWcbW vertex which enters the pomeron-exchange representation (b) for the generic 
Kancheli-Mueller discontinuity of the 4-body forward scattering amplitude is a highly 
nonlocal one, its properties will be discussed to more detail elsewhere. Upon integration 
over the momentum of the parton c one obtains the Kancheli-Mueller diagram (d) for the 
single-parton spectrum, for the nonlocal properties of the bWbW vertex in the diagram 
(d) see Sect .VI Dl. 
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V. NUCLEAR CASE: NONLINEAR fc ± FACTORIZATION 

The driving term of small-x nuclear structure functions and single-quark spectra from 
the 7* — > qq excitation off a nucleus were shown to take the familiar /c^-factorization 
form in terms of the collective nuclear Weizsacker- Williams gluon density 0(6, x A , k) per 
unit area in the impact parameter plane as defined in [2, 12, 13, 45]: 

(26) 

It satisfies the sum rule 

J d 2 K<p(b, x,k) = 1- exp[~(7 (x)T(b)] = 1 - S abs (b) . (27) 

where a (x) = J d 2 nf{x, n) is the dipole cross section for large dipoles. It is conveniently 
interpreted in the spirit of additive quark counting for soft hadronic processes as a (x) = 
o" 9 at(x) + u qN (x) = 2a q p N {x), so that the meaning of the factor 

S ahs (b) = exp[~a (x)T(b)] = exp[-af N (x)T(b)] (28) 

as representing intranuclear absorption becomes evident. 

Eqs. (17), (19) and (20) allow a straightforward calculation of the nuclear single-jet 
spectra. Nice analytic forms of the nuclear /c^-factorization are obtained in the large- N c 
limit, the 1/N% corrections can readily be derived following the considerations in [2] and 
will not be discussed here. The crucial point about the large- N c approximation is that 

a 3 (x, zr', r) =>- a qq (x, zr' - r) + cr q g(x, zr') , 

a gg (x, z(r - r')) => 2a qq -(x, z(r - r')) (29) 

so that 

T A [b,a 3 (zr',r)] = T A [b, a qq -(x, zr' - r)] + T A [b, <r qq -{x, zr')] 
- T A [b, a qq (x, zr' - r)} ■ T A [b, a qq (x, zr')} , 
T A [b,a gg (x,z(r-r'))] = 2T A [b, a qq (x, z(r - r'))] - T 2 A [b, a qq -(x, z(r - r'))] . (30) 



r A [b,a qq (x,r)\ = d K<j>(b,x,K 



exp[iH,r\ 



18 



FIG. 6: (Color online) Exchange by a nuclear pomeron IP^ in color dipole-nucleus scattering 
as a sum of multiple two-gluon pomeron exchanges. 

Then after some algebra we obtain 

(27i) 2 da A (g* -> QQ) _ 
dzd 2 pd 2 b 

Sabs(b) ■ J d 2 K(f)(b,x A ,K){\y(z,p) -^(z,p + zk)\ 2 + \V(z,p + k) - V(z,p + zk)\ 2 } 
+ J d 2 Kid 2 K20(6, xa, «i)0(b, %a, K 2 )\^(z,p + k 2 ) - ^(z,p + zk\ + zn 2 )\ 2 . (31) 

One can associate with (j)(b,XA, k>x) the exchange of a nuclear pomeron IP^ which as 
indicated in Fig. 6 sums multiple exchanges of the pomeron IP in interactions with 
nuclear targets ([2, 10], see also Sect. VI B). The first term is an exact counterpart 
of the free-nucleon result: in terms of the collective nuclear gluon density 0(6, Xa,k) 
it has exactly the same linear /c^-factorization form as the free nucleon result (24) in 
terms of f(xA,K>)- As such it corresponds to the Kancheli- Mueller diagram of Fig. 7a. 
Despite the similarity of diagrams of Fig. 7a and Fig. 5d, the two inclusive spectra are 
different because IP^ ^ IP. Furthermore, this contribution to the nuclear spectrum is 
suppressed by the absorption factor S a b s (b) and is entirely negligible for heavy nuclei. 
In this case only the second term survives, which is a manifestly nonlinear functional of 
the collective nuclear glue (j)(b,XA, k). It can be associated with the Kancheli-Mueller 
diagram of Fig. 7b. This completes our demonstration of the fc^-factorization breaking 
for single-jet spectra: the pattern of fcx-factorization changes from the conventional 
linear at the periphery of a nucleus to the nonlinear - quadratic - one for small impact 
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FIG. 7: (Color online) The Kancheli-Mueller diagrams for the mid-rapidity single-particle 
spectrum form gg subcollisions in pA collisions. The double-zigzag lines describe an exchange 
by the nuclear pomeron IP^- The open-charm spectrum in pp collisions is dominated by the 
diagram (a) with IP^ — > IP, see Fig. 5, the open charm spectrum in central pA collisions is 
dominated by the diagram (b) and the gluon jet production is dominated by the diagram (c). 

parameters. In retrospect, the finding of linear fc^-factorization for the single quark jet 
spectrum in 7* — > qq excitation can be traced back to the beam photon being colorless, 
by which the single-jet spectrum becomes an abelian observable [2]. In contrast to that, 
the origin of the nonlinear term in (31) is in the nonabelian nature of the 3-parton cross 
section (19), (29). 

The ^-dependence of the jet spectra enters through the wave-functions. To this end 
notice the strikingly different wave-function structure of the free-nucleon cross section 
(25) and its nuclear counterpart (31). This suggests a nontrivial dependence of the p- 
spectra on the longitudinal variable z which will be different for the free-nucleon and 
nuclear targets. Vice versa, the z-distribution will vary with p and the target mass 
number which can be called the LPM effect for open charm production. 
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VI. THE CRONIN EFFECT FOR OPEN CHARM 



A. The charm p-distribution: nuclear vs. free nucleon targets 

The impact of nonlinear /c^-factorization is quantified by the ratio of the free-nucleon 
and nuclear single-jet spectra, 

da A 



Ra/n(z;x a ,p) 



Ada 



N 



This ratio exhibits both shadowing, R A /n(z; x A ,p) < 1, and antishadowing, 
Ra/n(z; xa,p) > 1, properties, in the literature the latter is usually referred to as the 
Cronin effect. Our discussion of the Cronin effect is to a large extent based on the ideas 
and technique of the pQCD derivation of nuclear properties of the coherent diffractive 
breakup of pions into hard dijets [12, 13]. 

The principal points are best seen for central impact parameters, i.e., for S a bs(b) — > 0: 

(2n) 2 da A (g* -> QQ) _ 
dzd 2 pd 2 b 

J d 2 K 1 d 2 K 2 (j)(b,X A ,K 1 )(f)(b,X A ,K2)\^(z,p+K2)-^(z,p + ZK 1 +ZK2)\ 2 - (32) 

In all the comparisons of the free-nucleon and nuclear target observables we use the 
large- N c approximation without further notice. 

The results simplify greatly in two interesting cases. We start with the limit of slow 
quarks, z <C 1, when the free nucleon spectrum becomes 

2{27c) 2 da N (g* -> QQ) 



dzd 2 p 



Z<1 



J d 2 Kf(x A , K,)\y(z,p)-y(z,p+K)\ 2 . (33) 



In the limit of 2; C 1 the entire dependence of the nuclear spectrum (32) on Ki is 
confined to the normalization integral J d 2 Ki 0(6, x A , Ki) = 1 — S a b s (b) what lowers the 
nonlinearity of the nuclear cross section. Then, in conjunction with the linear term, 
oc S a b s (b), in (31), we find a result applicable at all impact parameters: 

(2n) 2 da A (g* - QQ)_ 

Z<1 



dzd 2 pd 2 b 



= J d 2 K<f>(b,x A ,K)\V(z,p + k) -*(z,p)\ 2 . (34) 
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We evidently recovered the linear /c_|_-factorization result (33), now in terms of the nu- 
clear collective glue 0(6, x a ,k). The physics behind this finding is the above discussed 
abelianization of g — > QQ at z <C 1, by which it becomes equivalent to 7* — > QQ. 

The second simple result is found for leading quarks with z — > 1 . Here the free-nucleon 
cross section again takes the form (33), whereas the nuclear single-jet spectrum can be 
cast in the convolution form 



da A (g* - QQ) 



dzd 2 pd 2 b 



Z->1 



y d 2 K(f)(b, x A , k) 



dcr A (x, p — k) 



+ S abs (b) 



da A (x,p) 



Z<1 



dzd 2 pd 2 b 



.(35) 



2<1 



<i,2<i 2 p<i 2 6 

The leading term in (35) is a manifestly nonlinear - quadratic - functional of 0(6, x A , k). 
Evidently, such a convolution would enhance the nuclear effects compared to the linear 
/c^-factorization formula (34). Arguably, when z varies from z — to 1, the p-dependence 
of the nuclear charm spectra shall interpolate between the ones given by equations (34) 
and (35). 



B. Properties of the collective nuclear gluon density 

The fundamental quantity of our formalism is the collective nuclear glue 0(6, xa,k) 
per unit area in the impact parameter plane [2, 12, 13]. It is characterized by a new 
scale - the saturation scale, 

47T 2 

Q\(b, x A ) w — a s (Q 2 A )G(x A , Q 2 A )T(b) , (36) 

a related quantity appears in other approaches to nuclear parton densities [35, 46]. At 
soft transverse momenta, n 2 < Q 2 A (b, xa), the nuclear collective glue exhibits a plateau 
which signals the saturation of parton densities. A useful parameterization in the plateau 
region is 

7T (k 2 + Q A (b,x A )) 2 

For hard gluons, k 2 >,Q 2 A (b,x A ), it is convenient to cite the unintegrated gluon density 
per bound nucleon: 

/a(6, xa-i k) = lm 
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1 + 



1 + 



o;5(K 2 )G(a; A , K 2 ) Q^(6,xa' 



z 

2 a5(gi)G(x A ,Q^) 
2ttV 



a 5 (^)G(:r A , K 2 )T(6) 



= /(xi, k)[1 + A ht (6,x^,k)] , 



(3* 



where 7 2 is an exponent of the large-/* 2 tail f(xA,K) oc (k 2 )" 7 . Here we show 
explicitly the leading nuclear higher twist correction, A#t(6, xa, k), which gives rise to 
a nuclear antishadowing effect [2, 12, 13]. To this end, the saturated gluon density (37) 
can be regarded as a resummation of all higher twist terms ~ (Q 2 A (b,x A ) / n 2 ) n . Note 
the very strong nuclear suppression of the gluon density per bound nucleon f A (b, xa, k) 
in the soft plateau region, n 2 <Q 2 A {b,x A ): 

2 1 



f A (b, x A , k) 



A -2/3 



(39) 



nQ\(b,x A )T(b) T 2 (b) 

The explicit expansion for (f)(b,x A ,K) in terms of the collective gluon density of j 
overlapping nucleons in the Lorentz-contracted nucleus [2, 12, 13] is 

<P(b,x A ,K) = S abs (b)j:j- (p^-J -f U) (^K) (40) 
where f^\x, n) is the j-fold convolution, 

fV\x, k) = f 6(£ «i - «) II fix, *i)d 2 Ki ■ (41) 
J 1=1 1=1 

Finally, we notice that (f)(b,x A , k) satisfies the important sum rule (see also [26, 33]) 



j2 2 



[b, x a ,k) - -T(b)f(x A ,K) 



0, 



(42) 



which readily follows from the property that F A [b,a(x,r)] — ^T(b)a(x,r) vanishes oc 



(T lx, v) as r 



0. The effect of nuclear suppression for nuclear soft gluons, see eqs. 



(37) and (39), is compensated for by nuclear antishadowing (38) for hard glue. 
The familiar integrated gluon density equals 



G(X,Q 2 ) = ^-J rfW /(*,«) 



(43) 
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its nuclear counterpart per bound nucleon, G A (b, x A , «^ 2 ), can be defined likewise in terms 
of f A (b,x A , k). Then the corollary of the sum rule (42) is that for Q 2 3> Q A (b,x A ) 

G A (b, x A , Q 2 ) - G(b, x A , Q 2 ) = - / d 2 K n 2 f(x A , n)A HT (b, x A , k) 

Jq 2 

2 1 2 n 2 a s (Q 2 ) a s (Q 2 )G(x A ,Q 2 ) Q 2 A (b,x A ) dG(x A ,Q 2 ) 



N c a s (Q 2 A (b,x A ))G(x A ,Q A (b,x A )) Q 2 dlogQ 2 

7 2 vr 2 a s (g 2 ) a s (Q 2 )Q 2 A(b, x A ) dG 2 (x A , Q 2 ) 



N c a s (Q 2 A (b,x A ))G(x A ,Q 2 A (b,x A )) dQ< 
where we made use of the expansion (38). 



(44) 



C. The Cronin effect: nuclear shadowing for p 2 <Q\(b,x A ) 

For the sake of illustration we start with the case of z 1 and assume that the satura- 
tion scale is so large that even for ttiq « p 2 « Q A (b, x A ) the jets are still perturbatively 
hard. Because the dominant contribution to (34) comes from k 2 ~ Q 2 A (b,x A ), the ex- 
pansion (B.6) of Appendix B simplifies further, and the k integration can be carried out 
explicitly 

(2n) 2 da A (g* - QQ) 

z«l V 

= [1 - S ahs (b)] ■ 2T F a s ■ — 2 . (45) 



dzd 2 pd 2 b 



2T F a s ■ — J d 2 K(f)(b,x A ,K) 



It is entirely dominated by the contribution from n 2 > p 2 , which in the language of 
evolution amounts to the anti-DGLAP splitting of collective nuclear gluons with large 
transverse momentum k into quarks with small transverse momentum p. 
The integration over all impact parameters gives 

2 J d 2 b[l - S abs (b)} = 2y>b{l - exp [-I (7o (x A )T(6)]} = a A [a (x A )\ , (46) 

which is the familiar Glauber-Gribov nuclear total cross section 04 [a ] for a projectile 
with the free-nucleon cross section a , and 

°Mf ~ A ~ m < 1 (47) 
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A marginal caveat is that for the fixed quark-jet momentum the inequality p 2 <C 
Q\{b,XA) does not hold uniformly over all impact parameters b, it breaks down for the 
most peripheral interactions, the peripheral contribution does not affect the shadowing 
property Ra/n(z; xa,p) < 1, though. 

Now turn to the opposite limit of z — > 1. Here the free-nucleon cross section is the 
same as for z <C 1. The convolution representation (35) for the nuclear spectrum is 
already suggestive of stronger nuclear effects. The analysis by the straightforward use of 
(32) in which one can neglect p in the arguments of the wave functions is still simpler: 

(27r) 2 da A (g* -> QQ) 

dzcPpcPb z _ yl 

2T F a s J dWMft x A , KxMb, x A , ^ {K 2 + m 2 Q){{ ^ +K2)2 + m 2 Q y (48) 

The integral in the r.h.s. will be dominated by the two logarithmic contributions from 
(i) k 2 , < k\ and (ii) (re 2 + ^i) 2 < n\: 

J d 2 Kl d 2 ^(b, XA , mMb, x A , ^) { ^ + ml){{ K 1 \ K2 y + ml) 

/r i 
d 2 K x (f)(b,x A , Ki)4>(b,x A ,0) + (f) (2) (b,x A , «i) log — J- 

« tt[0(6, a*, 0) + 0^ (b, x A , 0)] log (49) 



to. 



where <fi^(b, xa,k) is the convolution 

(2) (&,x a ,k) = (4>®(j)){b,XA,K). (50) 

For heavy nuclei (p^(b, xa,k>) has the same meaning as (/>(b,XA,K) but for twice larger 
T(b), what counts for the purposes of the present discussion is that 

[0(6, x A , 0) + 0( 2 ) (6, x A , 0)] « — -J- oc A' 1 / 3 (51) 

irQ A {b,x A ) 

The final result for the nuclear spectrum at z — > 1 is 

(27r) 2 <M<?* ^ QQ) 2T>a 5 (Q 2 A (6,x A )) Qjfo,^) 

10 S ~~2 



dzd 2 pd 2 b 



Q 2 A (b,x A ) * m 2 Q 
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The first dramatic effect is that the p-spectrum of fast quarks will be flat for p 2 < 
Q A (b,x A ) which must be contrasted to the oc 1/p 2 free-nucleon spectrum (45). The 
second dramatic effect is a strong enhancement of the shadowing effect. The result of 
the impact parameter integration can be evaluated as 

vr / d 2 b[^x A ^) +^\b,x A ,0)} « ^ XA \ (53) 

J \IJ A {0,X A )) 

where {Q 2 A (b,x A )) is an average saturation scale. This gives the shadowing effect 

d / i \ <7a[<to(x a )] P 2 , Q A (b,x A ) A - 2 /3 (KA \ 

Ra/n(z -> l;x A ,p) ~ — 7— ( — v- ■ jT^ajZ TTV ' lo S 2 K A ■ 54 

^0(2:4) ((Q A (b,XA))) m 2 Q 

The analytic formula for the transition from small to large z can be worked out making 

use of the explicit parameterization (36). The gross features of this transition are well 

reproduced by the interpolation formula 



c / \ cr A [a (x A )} 

R A/N (z;x A ,p)K a<Xa) 

applicable for p 2 <(Q A {b,x A )). 



P 2 , ]n (Q 2 A (b,x A )) 



p 2 + z 2 (Q A (b,x A )) (Q 2 A (b,x A )) ta m, 



2 



(55) 



D. The Cronin effect: hard quark-jets, p 2 >^Q\(b,x A ) 
1. Open charm from the direct and resolved gluon interactions 

The above /cj_-factorization formulas are exact and can be directly applied to the 
calculation of the jet spectra. In order to make a contact with the more familiar treatment 
of photoproduction of jets, here we present the LLp 2 decomposition of the jet spectra into 
the direct and resolved interactions of the incident parton a. Such a decomposition will 
make more obvious the origin of the Cronin effect, although in the practical calculations 
there is no need to resort to the LLp 2 approximation. 

For hard jets we can use the large-p approximation (B.6) of Appendix B. Then the 
hard single-jet spectrum for the free-nucleon target takes the form 

da N (g* -> QQ) _ 
dzd 2 p 

26 



T F a s (p 2 )[z 2 + (1 - zf 



J d 2 KK 2 f(x A ,K) 



(2n) 2 

whereas for the heavy nuclear target 

da A (g* - QQ) 2T F a s (p 2 ) [z 2 + (1 - zf] 



z 2 (1-z) 2 



p 2 (p + k,) 2 



(p + zk)'- 



, (56) 



dzd 2 pd 2 b (2tt) 2 

r l2 l2 .,. s (^Ki — (1 — z)k 2 ) 2 

X / d Kid K,2(p{b,XA,Kl)nb,XA:K 2 )-, ; ^7 ; ; r 

J [p + K 2 y(p + ZK! + ZK 2 ) 



2 



(57) 



For the sake of brevity here we suppressed the contribution from the linear term oc S a t> s (b) 
in (31). 

In the spirit of LLp 2 one can disentangle two distinct sources of the transverse mo- 
mentum of the quark-jet. In the first case the pQCD subprocess can be viewed as the 
large-p direct (unresolved) boson-gluon fusion g*g — > QQ where the exchanged gluon g 
has a small transverse momentum \k\ < \p\, the large p flows in the t-channel of the 
pQCD subprocess g*g — * QQ. Such a contribution can be evaluated as 



d<T N (g* -> QQ) 



dzd 2 p 



dir 



T F a s (p 2 )[z 2 + (1 - z) 2 ] 2 An 2 a s (p 2 ) 2 , 



(p 2 ) 2 N c 
2n[z 2 + (1 - z) 2 } A7r 2 a s (p 2 ) dQ(z,p) 
p 2 N c dp 2 



G(x A ,p 2 ) 



G(x A ,p 2 ). (58) 



First, one recovers the familiar collinear factorization proportionality to the target gluon 
density G(xa,P 2 )- Second, one readily identifies the unintegrated density of large-p 
quarks in the gluon, 

dQ(z,p) _ T F a s (p 2 )[z 2 + (l-z) 2 ] 

dp 2 2np 2 ■ 1 j 

One can say that the direct process probes the transverse momentum distribution of 

quarks in the beam gluon g*. 

The second contribution corresponds to the t-channel quark being close to the mass 

shell, (p + k) 2 ~ m 2 Q. In this case the large transverse momentum p of the quark-jet 

comes from the large transverse momentum of the exchanged gluon, k —p, in close 
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similarity to the pomeron splitting mechanism for production of hard diffractive dijets 
[47]. It also can be viewed as a jet from the resolved gluon interactions. All the target 
information is encoded in the unintegrated f(x A , «) for a free nucleon or <f>(b,XA,p) for 
a nucleus, and the process itself looks as the fusion Qg — > Q'. One must be careful with 
the fusion reinterpretation as the resolved contribution is only a part pf the more generic 
9*9 ~^ QQ an d singling it out only makes a unique sense to the LLp 2 approximation. 

The case of finite < z < 1 is a bit more involved and will be considered to more 
detail elsewhere. The two poles in (56), at (p + zn) 2 = and (p + n) 2 = 0, are well 
separated. At z — > 1 the two poles merge but the residue vanishes oc (1 — z) 2 . Here-below 
for the sake of illustration we focus on slow quark production, z <C 1. For the massive 
quarks upon the azimuthal averaging 

1 11 

(p - K) 2 + m 2 Q ) ^ ^(p2 _ K 2 _ m 2 Q y + Ap 2 m 2 Q =^ |p2 _ K 2| + $2 ' ( 60 ) 

where the infrared regulator of the dn 2 integration can be taken as 5 4 w p 2 m 2 Q . The 
corresponding contribution to the quark-jet cross section for the free-nucleon target can 
be evaluated as 



da N (g* -> QQ) 



dzd?p 



T F Tia s {p 2 ) . . p 2 
-^2-f^P)^ Q 

2tt 4n 2 a s {p 2 ) dG{x A ,p 2 ) 



Q(z«l,p). (61) 



p 2 A" c rfp 2 

The emerging logarithm times the splitting function can be identified with the integrated 
quark structure function Q(z l,p 2 ) of the beam gluon g* probed by the exchanged 
gluon g of virtuality Q 2 = p 2 . The resolved gluon mechanism probes the transverse 
momentum distribution of gluons in the target. 

In the more general case one must take the quark density Q(xb,P 2 ) where xb is the 
fraction of the lightcone momentum of the beam carried by the quark Q. Combining 
together (58) and (61) we find the spectrum of hard quark-jets which is proportional to 

(dQ(x B ,p 2 ) 2 2 dG(x A ,p 2 )\ d 2 2 

[ dp~ 2 G ( XA >P ) + Q( x b,P ) ^ J = -jjp[Q{x B ,P )G(x A ,P )] , (62) 
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which restores the beam-target symmetry. A similar observation has been made in [26]. 
The last form shows that in the general case the both mechanisms are of comparable 
importance. It also demonstrates the nonlocality properties of the 6IP6IP vertex in the 
Kancheli-Mueller diagram of Fig. 5d. The derivation of the nuclear counterparts of (58) 
and (61) for z <C 1 is straightforward and need not be repeated here. We proceed directly 
to the derivation of the Cronin effect. 



2. Antishadowing property of collective nuclear gluon density as an origin of the Cronin 
effect 

For a more accurate isolation of the Cronin effect for hard jets consider the nuclear 
excess cross section 

dAa A da A 1 

-tfb = M ~ 2 T{b)daN ■ (63) 
Again we start with z < 1. Suppressing common factors and making use of (60) we 
need to evaluate 

dAdA r „ k, 2 r , 1 , 



'A 

OC 



/ d " K \ 2_ n K? \ +5 ^ XA ^ ) - 2 7 W( a *'")] (64) 



dzd 2 pdb J \p 2 

A priori the sign of this quantity is not obvious because the function A0^(6, xa, «) = 
(f)(b,x A ,K<) — |T(6)/(x J 4, k) is negative valued for n 2 < Q 2 A (x A ,b) and positive valued 
in the antishadowing region described by eq. (57). We decompose n 2 /(\p 2 — n 2 \ + 5 2 ) as 
follows: 

9(p 2 - k 2 ) ■ % + 9(k 2 - p 2 ) + K(p 2 , k 2 ) . (65) 



\p 2 — K? | + 5 2 " p 2 

Here the first term isolates the direct boson-fusion contribution, whereas the collinear 
quark-pole contribution is contained entirely in K(p 2 , n 2 ) which can readily be shown to 
be a positive- valued function which vanishes ~ k 4 for n 2 <^ p 2 and ~ 1/n 2 for k 2 3> p 2 . 

The shadowing effect from the negative- valued A0^(6, x A , k) at small p 2 < Q\(b, xa) 
is concentrated in the contribution from the first term in (65) which by virtue of the sum 
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rule (41) can be represented as 
1 rP 2 



1 rP 1 f°° 

— dK 2 K 2 A(j) A (b,X A ,K) = -/ dK 2 K 2 A(j) A (b,X A ,K) 

p 2 Jo p 2 Jp 2 



(66) 



Consequently, the nuclear excess cross section is entirely calculable in terms of gluon 
densities in the hard region of n 2 > p 2 : 



2 dAa A 
p - — - — — oc 



1 r°° 

+ / dK 2 A(j) A (b, x a ,k) 
Jp 2 

+ - f d 2 KK(p 2 ,K 2 )A(p A ( y b,x A ,K) 

7T J 



(67) 



Now notice that in this hard region A<fi A (b, x A , n) oc 1/(k 2 ) 3 , see eq. (38), and all 
integrals in (67) are well converging ones. The collinear quark-pole contribution can 
be evaluated as it was done in (61), factoring out p 2 A(j) A (b,p), so that to logarithmic 
accuracy 



P JTynnz K A<f) A {b,p) 



■ 1 , P 2 

-l + ^ + log^ 



(68) 



dzd 2 pdb 

Following the derivation of eq. (61) we can associate log(p 2 /mg) with Q(z l,p) and 
to LLp 2 accuracy 

dAa A daw 



dzd 2 pdb dzd?p 
daw 



f A (b,x A ,p 2 ) _ ' 
f{xA,P 2 ) 



dzd?p 



T{b) 

T(b)A HT (b,x A ,p). 



(69) 



The result for the resolved interactions does not depend on the projectile parton a, 



[RA/N(Z < 1;ZA,P) - l]res = A HT (b,X A ,p) 



(70) 



while the dilution of the Cronin effect by the contribution from direct interactions to the 
free-nucleon is process dependent: 

d<7N(res) 



Ra/n(z < 1;xa,p) - 1 = 



dcr N (res) + da N (dir) 



[R A/N (z < l;x A ,p) ~ l]re«- (71) 
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Hereafter we focus on the projectile independent antishadowing for resolved interaction 
and suppress the subscript 'res'. For heavy nuclei 



and we obtain our final estimate for the antishadowing Cronin effect at large-p 2 



Here we adopted the same subscript Qg as in the familiar splitting functions. 

First, of the two components the direct g*g-fusion would have given nuclear shad- 
owing, see the first line in eq. (67), the antishadowing nuclear excess cross section is 
a feature of the resolved gluon interactions. Second, the fundamental point is that for 
hard jets the nuclear excess cross section follows directly from the nuclear antishadowing 
component A#t(&, x A ,p) of the collective nuclear gluon density [12, 13]. Third, the an- 
tishadowing rises with the nuclear mass number, R A / N (z; x A ,p) — 1 oc A 1 ^. Fourth, the 
antishadowing nuclear excess is a special quadratic functional of the gluon density in the 
proton - a product of the integrated and unintegrated gluon densities. Finally, for hard 
jets the antishadowing Cronin effect vanishes oc l/p 2 . Whereas the antishadowing for 
hard jets, Ra/n(z] x a , p) > 1, has been discussed by many authors under model assump- 
tions and specific parameterizations for the dipole cross section/unintegrated glue (e.g. 
[24-26, 33] and references therein), the above reported model independent derivation, 
and identification of the source, of antishadowing Cronin effect are new results. 

3. Variations of the Cronin effect from slow to leading charm 

The convolution representation (35) makes it obvious that the above findings on the 
Cronin effect are directly applicable to the opposite limiting case of z — > 1 as well. We 
recall that the free-nucleon spectra for z — > 1 and z C 1 are identical. The smearing of 
the decreasing spectrum (34) does obviously enhance the antishadowing effect at large p 
which can readily be seen as follows. The p-dependence of the spectra (58) and (61) is 




(72) 



A Qg (z < l;x A ,p) = R A/N (z < l;x A ,p) - 1 



9vr 7 2 a s (p 2 )-G(x A ,p 2 ) A 



(73) 



l.V r ' ' ~' p 2 R\ 
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driven by the factor 1/ (p 2 ) 2 and the enhancement of the Cronin effect will for the most 
part be due to the smearing of l/[(p — k) 2 ] 2 in (35). The large-p expansion 



1 + 4— (74} 
[(p-«) 2 ] 2 (p 2 ) 2 [ p 2 \ ' { V 

in which the azimuthal averaging is understood, gives an enhancement of Ra/n( z ^ 
1; xa,p) by the extra factor 

PA/N(P) = J d 2 K(f)(b,XA,K) 

a s (p 2 )G(x A ,p 2 ) Q 2 A (b,x A 



1 + 4^ 



1 + 2 



as(Q 2 A )G(xA,Q A ) 



l + A HT (b,x A ,p), (75) 



cf. eq. (38) for 7^2. Consequently, at z — > 1 the antishadowing Cronin effect will be 
twice stronger than that at z <C 1, 

Aq 9 (,2 -> 1;za,p) = Ra/n(z < 1;^a,p) • Pa/n(p) - 1 ~ 2A Qs (2; < 1;za,p) , (76) 

what nicely correlates with the stronger nuclear shadowing at small p. The same corre- 
lation is obvious from the sum rule which relates the p-integrated nuclear cross sections 
for z — > 1 and z <^ 1, 



da A (g* -> QQ) 



dzd 2 b 



da A (g* -> QQ) 



Z->1 



dzd 2 b 



(77) 



Z<1 



which readily follows from the convolution representation (35). 



^. Numerical estimates for the Cronin effect 

The antishadowing properties of the collective nuclear gluon density have been studied 
[12, 13] in connection with the A-dependence of coherent diffractive dijet production. 
The starting point was the unintegrated gluon density in the proton, J(x,k), adjusted 
to reproduce the 7*iV interactions from real photoproduction, Q 2 = 0, to DIS at large 
Q 2 [14]. The quantity of interest is the antishadowing higher twist effect in the collective 
nuclear gluon density for j overlapping nucleons which is given by the j-fold convolution 
f^\xA, k). This density per overlapping nucleon was represented in ref. [12] as 

fA(b,x A ,K) = J r*** = f(x A , k) [1 + (j - 1)A,(^, n)] . (78) 
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With reference to the Poissonian form of the multiple scattering expansion (40), in 
Figs. 8,9 we show Ac r0 nin( x A, k) — j • Aj(xa, k) for two values of xa- 

These results for Ac r <min{%A, «) can be used to evaluate A(xa,p) = Ra/n(xa,p) — 1 
for different jet production subprocesses and for different nuclei: 

A Qg (z < 1; x A , p) ~ A Cronin (x A , k) (79) 

For the 196 Pt nucleus, which is not any different from gold - the favorite of RHIC, the 
average value of j in the expansion (40) is (j) ~ 4. At xa = 0.01 the peak value of the 
antishadowing effect, Aq 9 (z <C 1;xa,p) ~ 0.7, will be reached at p ~ (2 -i- 2.5) GeV. 
It is important that although the numerical studies in [2] suggest not so large Q\ < 1 
GeV 2 , the peak of the antishadowing effect is placed at transverse momenta p which 
are sufficiently large to justify the pQCD evaluations. For quark jets with z — > 1 the 
antishadowing effect will be twice stronger. The relationship between Ac r <min(^A, t) and 
A gg (z; xa,p), A gq (z; xa,p) for gluon jets will be discussed in Sections VII B and VIII B, 
respectively. 

The x^-dependence of the Cronin effect is noteworthy. At xa = 0.001 the peak value of 
the antishadowing effect, Aq 9 (z <C 1; xa,p) ~ 0.25, will be reached at p ~ (3^-3.5) GeV. 
This x-dependence of the position and magnitude of the peak value of the antishadowing 
effect is driven by three effects. The first one is the x-dependence of Q A { x i b) which is 
similar to that of the real photoabsorption cross section and is rather weak, ~ x~ 0A . The 
second source is the variation of the /^-dependence of f(x, k) as x decreases. Indeed, the 
phenomenological studies of ref. [14] show that while the exponent 7 ~ 2 is appropriate 
for xa = 0.01, at xa = 0.001 the /^-dependence of the unintegrated gluon density in the 
region of k 2 of interest corresponds to 7 ~ 1.7 which lowers substantially the magnitude, 
and shifts to larger n 2 the onset of the antishadowing regime (38). Finally, because of 
the shift to larger p 2 with decreasing x the peak value of A Cr onin(p) is suppressed by 
the 1/p 2 dependence of the antishadowing contribution, see eq. (73). This explains the 
numerical findings of Ref. [12]. 

We recall that the shown numerical results are only for the resolved interactions, any 
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p [GeV] 

FIG. 8: (Color online) The antishadowing Cronin effect in jet production off a heavy nucleus at 
x = 0.01 as a function of the transverse momentum p of the jet parton for several multiplicities 
of overlapping nucleons j in the Lorentz-contracted nucleus. The production of slow quarks 
off nucleus with mass number A = 200 corresponds to j = 4; the production of slow gluons off 
nucleus with mass number A = 200 corresponds to j = 8. 

comparison with the experimental data requires an application of the dilution factor (71). 
Also, the momentum of the observed hadrons is only a fraction of the jet momentum 
what places the Cronin peak for hadrons at smaller values of p than for jets. 
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p [GeV] 



FIG. 9: (Color online) The antishadowing Cronin effect in jet production off a heavy nucleus at 
x = 0.001 as a function of the transverse momentum p of the jet parton for several multiplicities 
of overlapping nucleons j in the Lorentz-contracted nucleus. The production of slow quarks 
off nucleus with mass number A = 200 corresponds to j = 4; the production of slow gluons off 
nucleus with mass number A = 200 corresponds to j = 8. 

VII. NONLINEAR /cx-FACTORIZATION FOR PRODUCTION OF GLUONS 
OFF NUCLEI 

A. The gluon spectra for the free-nucleon and nuclear targets 

The excitation of the gluon g — > gg offers a still more remarkable example of breaking 
of linear /^-factorization. In this case a = b = c = g and we take the shorthand notations 
z = z g and p = p . The 3-parton cross section for generic three-gluon state equals 




C A 



[<j q q{x, p) + a q g(x, p - r) + a q g(x, r)) 




2C F 
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The dipole cross section structure in eq. (16) takes the form 

C 

+ cr qq -(x, r') + cr qq (x, zr) + a qq (x, r' - zv) 

- 2<j gq -(x, r-r')- 2a g g(x, z(r - r'))} (81) 

Repeating the analysis of sections III and IV, for the free-nucleon case one readily finds 
the /c_i_-factorization representation for the gluon-jet cross section 

+ \V(z,p+ k) - ^(z,p + zk)\ 2 + \^{z,p) - ^(z,p + k)| 2 } . (82) 

In agreement with the generic arguments of section IV, the abelian limit is recovered for 
z < 1: 



dcr N (g* -> gg) 



dzd 2 p 



2<1 



2(2tt) 2 



J d 2 Kf(x A , k) ^ • \*(z, p) - *(z, p + k) I 2 . (83) 



i / ... c A 

The same abelian limit is recovered also for z — > 1, but this equality of two limiting cases 
will again be upset for nuclear targets. 

Following the discussion in section VI D, one can readily decompose the spectrum 
(83) into the contributions from the direct and resolved gluon interactions, we shall skip 
those details here. 

The excitation of gluons into gluons, g — > gg, is an exceptional case for which exact 
formulas for a nuclear target derive without invoking the large- N c approximation. Our 
result for the nuclear gluon jet spectrum reads 

w ^ar w) -[^w]'-/^^-») 

x{|tf(2,p) - V(z,p + zk)\ 2 + \V(z,p + K) -^{z,p + zk)\ 2 + \V(z,p) -V(z,p + k)| 2 } 

+ S abl( b ) J d 2 K 1 d 2 K 2 (f)g(b, X A ,Ki)(() g (b,X A ,K2) 

x[\m( Z ,p+ m) - ^(Z,p + ZK 2 )\ 2 + \V(Z,P+ Ki + K 2 ) - ®(Z,P+ ZK 2 )\ 2 
+ \-V(z, P + «i) - V( Z ,p + z(Ki + K 2 ))| 2 } 

+ J d 2 K 1 d 2 K 2 d 2 K 3 (j)g(b, XA,K 1 )^ g (b,XA,K 2 )(l>g{b,XA,K3) 

x|tf(z,p + Ki + K 3 ) -^(z,p + z(k 2 + k 3 ))\ 2 . (84) 
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Here 4> g {b,x A , «) is the collective nuclear glue which now must be evaluated from 

T A [b,^-<J q q(x A ,r)} = J d 2 K(f) g (b,XA, k) 1 - exp[inr] , 

S ( £ s (b) = exp[-i ■ ^- ■ a (x A )T(b)} . (85) 

in which the nonabelian factor Ca/2Cf enters manifestly. The emergence of a new 
collective gluon density (p g (b,XA, k) 7^ <f>(b,XA,K) illustrates nicely the point made in 
[2] that the collective nuclear glue must be described by a density matrix in color space 
rather than by a universal scalar function. Correspondingly, the nuclear pomeron W 9 A 
associated with <p g (b, xa, k) is different from IP^ defined for charm production off nuclei 
in Sect. V. The saturation scale for the gluon density 4> g (b, xa, k) differs from Q\{b, xa) 
by the factor Ca/ICf- Also, for hard gluons 



f A 9) (b,x A ,K) = f(x A ,K) - 



1 + 7^r^HT(b, x a ,k) 



(86) 



2C F 

In close similarity to open charm production, the first term in eq.(84) which satisfies 
linear /cj_-factorization (82) subject to the substitution of the free-nucleon f(xA, k) by 
the collective nuclear glue (f> g (b, xa, k), is suppressed by the square of the nuclear absorp- 
tion factor, S$ s (b) . The second term, which is quadratic in the nuclear gluon density 
(f) g (b,XA, k), is suppressed by s£ s (b). What survives for strongly absorbing nuclei is 
the third term which is a cubic functional of collective nuclear glue <p g (b,XA, k) and the 
pattern of /cx-factorization for single gluon jet production changes from the conventional 
linear one for peripheral impact parameters to the cubic one for small impact param- 
eters. The three components of the gluon spectrum can be associated with the three 
Kancheli-Muller diagrams of Fig. 7. A comparison of (82) and (84) shows that the nu- 
clear cross section cannot be represented as an expansion in multiple convolutions of the 
free nucleon cross section as was suggested in [28]. The z-dependence of the integrands 
of the quadratic and cubic terms is manifestly different from that of the free-nucleon 
cross section. The difference between the the free-nucleon and nuclear gluon densities 
entails different ^-distributions for free-nucleon and nuclear targets even in the linear 



37 



term. Consequently, all the three components of the nuclear spectrum will exhibit the 
p-dependent LPM effect. 

The cubic nonlinearity in terms of the collective gluon density is a pertinent feature 
of the radiation of gluons with finite z. Some interesting simplifications are found in the 
limiting case of z <C 1. In this case in the quadratic and cubic terms one will encounter 
the quantities of the from 

J d 2 K 1 d 2 K 2 <frg(b,X A ,K 1 )(j) g (b,X A ,K2)\^(z,p + K, 1 + K 2 ) ~ ^(z,p)\ 2 

= J d 2 ^f (b,x A ,K)\V(z, P + K)-V(z,p)\ 2 . (87) 



In terms of the convolution 4> g 2 \b, x A , k) the resulting single jet spectrum can be cast in 
the form 

(27r) 2 da A (g* -> gg) 



C f 

T A b , yrVqqixi r)] = / d 2 K(f) gg (b, x, k) \l - exp[inr]] (90) 



dzd 2 pd 2 b 

J d 2 K[2S%l(b)<f> g (b,x A ,K)+0 (b,x A ,K)} .\*(z,p + K,)-*(z,p)\ 2 , (88) 

and the triple- IP^ exchange Kancheli- Mueller diagram of Fig. 7c simplifies to the double- 
IP^ exchange diagram of Fig. 7b. Finally, we observe that 

<P gg (b,x A ,K)=2S { ±(b)<P g (b,x A ,K) + <pf\b,x A ,K,) (89) 

is precisely still another collective nuclear gluon density defined in terms of the gluon- 
gluon dipole cross section: 

(_ 

To this end, eq. (88) is analogous to the linear /cj_-factorization for the free nucleon 
case, see (83), but the free nucleon glue f(x A ,n) is substituted by 4> gg (b, x A , k), which 
is still a nonlinear -quadratic- functional of the collective nuclear glue (f> g (b, x A , k), 
and with which one must associate still another nuclear pomeron JP g A g . In terms of 
(p gg {b,x A , k) the transverse momentum spectrum of gluon jets with 2; < 1 satisfies the 
linear /cj_-factorization and is described by the Kancheli-Mueller diagram of Fig. 7a with 
IP — > IP^f . Our final result for z 1 coincides in its color dipole form with the one 
cited in [20], see also a discussion in [2], all the results for finite z are new. 
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B. The Cronin effect for gluons 



The implications of (90) for the Cronin effect for gluon-jets are straightforward. The 
principal change is that for soft gluons the collective nuclear gluon density (f> gg (b,XA, «0 
is obtained from <f>(b,XA, «) by the substitution 

Q 2 A (b, x A ) Q% gg (b, x A ) » ^Q\(b, x A ) , (91) 

so that the plateau in (f) gg (b,XA, n) will be broader than in <f)(b, xa, k) by the factor 
Ca/C f . For hard gluons 



f A "\b, x A , k) = f(x A , k) ■ ^ 



C 'a 

1 + —A HT (b, x a ,k 2 ) 



(92) 

Nuclear shadowing for p 2 < Q\ gg (b, x A ) will be stronger than for the quark jets: 



R% N {z;x A ,p) = 2^ a) ~ ^R% N (z;x A ,p) (93) 

Evidently the antishadowing effect A gg (z; x A , p) > at p 2 > Q\ gg {b, x A ) will persist for 
gluon-jets too. Following the analysis of section VI D, one would readily find that the 
nuclear antishadowing effect will be dominated by the resolved gluon interactions, and 

C 

A gg (z;x A ,p) w -^-A Qg (z;x A ,p) . (94) 

U F 

In view of the larger saturation scale (91) the position of the peak value of the anti- 
shadowing effect for gluon jets, A gg (z; x A ,p), will be placed at higher p 2 than that for 
quark jets. For semiquantitative estimates we again invoke the numerical results for 
Acronin{x A) p) shown in Fig. 4 of Ref. [12]. For heavy nuclei like gold or platinum 
(j) g ps §^0) ~ 9. At x A = 0.01 quite a large peak value A gg (z;x A ,p) ps 1.3 will 
be reached at p ~ 3 GeV. At x A = 0.001 the peak value A gg (z; x A ,p) ~ 0.45 will be 
reached at p ~ 5 GeV. Evidently, the fragmentation of jets will shift the peak position 
for observed hadrons to lower values of p. The antishadowing properties of the spectrum 
of leading gluons with z — > 1 are readily derived following the discussion of quark jets 
in section VI D. The full treatment of virtual corrections to the production of leading 
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gluons from the lower order pQCD processes is beyond the scope of the present commu- 
nication, here we only notice that they are arguably small if the transverse momentum 
p of the observed leading gluon is much larger than the transverse momentum of the 
beam-gluon g*. The leading gluon version of eq. (84) demonstrates how the pertinent 
nonlinearity at z — > 1 changes the antishadowing Cronin effect compared to the above 
discussed case of z <C 1. Define first the two special nuclear spectra 
(2n) 2 da { f 



dzd 2 pd 2 b 
(2n) 2 da ( f ) 
dzd 2 pd 2 b 



J d 2 K(j) g (b,X A ,K)\y(z,p)-y(z,p+K)\ 2 (95) 

J d 2 n l d 2 n 2 (t)g{b,x A ,n 1 )(t) g {b, x A , k 2 ) \V(z, p + «i) - V(z,p + k 2 )\ 2 , 

(96) 



in terms of which 

dcr A (g* -> gg) 



dzd?pd 2 b 



Z->1 

, , lal ) \ , ■ . iti2) \ . ^, 



2 da { l ] (g* -> gg) (da { f ] <g> <t> g ) (p) {da { f ] ® 4> g ) I 

2Sabs{b) dzd 2 pd 2 b + 2Sabs{b) dzd 2 pd 2 b + dzd 2 pd 2 b ■ (97) 
The small and large-p properties of the first two terms in (97) have already been studied 
in section VI. After averaging over the azimuthal angles of n lt2 : 

*(2!,P+Kl) -^(Z,p + K 2 )\ 2 

1 \ / 1 \ 2 



(p+Kl) 2/ X (P + K 2 ) 2 ' p 



+ { ¥ l—^)-^P 2 -^P 2 -^ 

1 2 



» 2 -(Hp 2 *)0(p 2 h.z). (98) 
\p 2 - k(\ + o 2 p 2 

where we made a provision for the K\ <-> /«2 symmetry of the integrand. At small p 
the last term in (98) can safely be neglected and the small-p behaviour of da^ will be 
similar to that of the spectrum of leading quarks (52) subject to swapping the g — > QQ 
splitting function for the g — > gg one as described in Appendix B. For the derivation of 
large-p behaviour a more convenient representation is 

^p + ^-^P + ^l 2 )) = 7^(p 2 - >A) + - p 2 ) 



- K(p 2 , nj) + - 2 9{p 2 - ni)6(4 ~ P 2 ) 
p 
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(99) 



The logarithmic pole contribution is absorbed in K(p 2 ,kI) which has the same prop- 
erties as K(p 2 ,k\) in the expansion (65). Then the discussion of nuclear shadow- 
ing/antishadowing properties of these three contributions will be entirely identical to 
that for the quark spectrum in section VI D. Following the common wisdom the loga- 
rithm log(p 2 / ixq) in the resolved gluon contribution, where \iq is the appropriate infrared 
cutoff for gluons, must be associated with the integrated density of gluons in the beam 
gluon, G B (z,p 2 ). The new item in the expansion (99) is the last term; one can readily 
verify that its contribution is short of the logp 2 and can be neglected compared to that 
from K(p 2 , k 2 ). Then, to the leading logp 2 approximation, 



dzd 2 pdb 



dcr N (g* -> gg) 



dzd 2 p 
da N (g* -> gg) 



dzd 2 p 



T{b) 



T{b) 



2C F f ( A 9 \b,x A ,p 2 ) 
C A }'{xa,P 2 ) 
C A 



1 + 



2C, 



-A HT (b, x A ,p) 



(100) 



The dominant contribution to (97) comes from the last, nonshadowed, term in which the 
spectrum (100) is convoluted with <f) g (b,x A , «)■ The convolution effect has been derived 
in section VI D, see eq. (69), the generalization to our case gives the extra factor 

C A 



1 + 



2C, 



-A HT (b : x A: p) 



(101) 



and 



dAa A (g* -> gg) 



dzd 2 pd 2 b 



dcr N (g* -> gg) 



T(b)-^A HT (b,x A ,p). (102) 



z~>i dzd 2 p 

A remarkable finding is that for this component of the gluon spectrum the antishadowing 
Cronin effect does not change from z < 1 to z -> 1. One can attribute that to the 
observation that although the nonlinear nuclear /c^-factorization in the two limits is of 
quite distinct form, the nonlinearity in the both cases is the same - a quadratic one. 
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VIII. GLUON RADIATION OFF QUARKS AND THE SPECTRUM OF 
QUARKS FROM q* -► gq 



A. Linear factorization for gluon radiation off free nucleons 

For completeness we show here the spectrum of gluons radiated from a fast quark 
that propagated the nucleus, as well as the spectrum of quarks after having shaken off 
their radiation cloud. 

Here a = q*(b), b = g(b g ), c = q(b q ), where we also indicated the corresponding 
impact parameters. The lightcone momentum fraction carried by the gluon is z g . Because 
the transverse momentum distributions of produced gluons and scattered quarks are 
different, for the sake of clarity we done them by p g and p , respectively. 

The relevant dipole distances are: 

b g -b q = r,b'-b q = z g r',b g -b' g = r-r',b-b' = z g (r- r') . (103) 

Again we start from the free nucleon case. The differential cross section is written as 
2{2Ti) 2 da N {q* -> g(p g )q) 



dZgffp 



= J d 2 rd 2 r' exp[ip g (r - r')]^(z g , r)^*(z g , r') 

x {<7 3 (r, z g r') + <r 3 (r', z g r) - a gg (r - r') - a qq -(z g (r - r'))} , 

(104) 



and, again, now in momentum space: 

2(2n) 2 da N (q* -> g{p g )q) _ 

dz g d 2 p g 
r C 

J d 2 Kf(x A , K ){^-[\^( Zg ,p g ) - V( Zg ,p g + K,)\ 2 + \V(z g , Pg + K ) - V(z g , Pg + Z g K)\ 2 
-\V(z g , Pg ) - y(z g , Pg + Z g K)\ 2 ] + \*(z g , Pg ) ~ ^(Z g ,p g + Zg K)\ 2 } . (105) 

Which of course is nothing else but the QCD-'Bethe-Heitler'-Bremsstrahlung-spectrum 
[8] albeit in a not so familiar representation in terms of lightcone wavefunctions. 
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B. Nuclear case: nonlinear k±_ factorization for gluon radiation 



We quote our result in the momentum space form: 

(2n) 2 da A (q* -> g(p g )q) _ 
dz g d 2 p g d 2 b 

Sabs{b) J d 2 K(f)(b,X A ,K){\V(z g ,p g ) -q{z g , Pg + K)\ 2 + \V{z g , Pg + K.) - V{z g , Pg + Z g K,)\ 
+ J d 2 K 1 d 2 K 2 <l)(b,X A ,H, 1 )(f)(b,X A ,K2)\^(z g , Pg + Z g K t ) ~^(z g , Pg + Ki + K 2 )\ 2 ■ (10 

Note the by now familiar decomposition into the absorbed linear and nonlinear - 
quadratic -terms. 

Evidently, the spectrum of soft gluons does not depend on the spin of the beam 
parton. Indeed, for soft gluons, z g <C 1, the result (106) simplifies to 

{2ix) 2 da A {q* -> g[p g )q) 

dz g d 2 Pg d 2 b zg<<1 

J d 2 K{2S abs {b)<f)(b,x A ,K) +<f) (2 \b, x A ,K)}\V{z g ,p g ) -V(z g , Pg + K)\ 2 

= J d 2 K(f) gg (b,x A ,K)\y(z g ,p g ) -V(z g ,p g + k)\ 2 , (107) 

which is a precise counterpart of eq. (88) for the g* — > gg breakup. In the large-iV c 
approximation 4> g (b, x A ,n) = (f)(b,x A ,n) and (j) <8> (f)j (b, x A , k) = 4> gg (b, x A , k). 

Consequently, the antishadowing Cronin effect for slow gluons from q* — > qg will be 
identical to that from g* — > gg, 

A gq( z g < l;x A ,Pg) = A gg(z g < l;x A , Pg ) = 2A Qg (z Q < V,x A ,P g ) ■ (108) 

For the numerical estimates for pAu collisions one must take Ac r onin(x A , Pg ) for j — 8 
in Figs. 8, 9. On the other hand a comparison of eqs. (106) and (31) shows that the 
spectrum of leading gluons with z g — > 1 from q* — > qg will have precisely the same 
convolution form (35) as the spectrum of leading quarks from g* — > QQ which in view 
of the findings of sections VI D and VII B entails 

A gq {z g -> l;x A , Pg ) = A Qg (z Q -> l;x A , Pg ) = A gg (z g -> l;x A , Pg ). (109) 

This finding suggests that the equality of antishadowing Cronin effect for leading and 
soft gluons does not depend on the beam parton be it a quark or a gluon. 
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C. p q spectrum of quarks from the breakup q* — > gq on a free nucleon 

We now turn to our last application, the spectrum of quarks from the excitation of 
the gg-Fock state of a fast quark. In the practically interesting case of forward quark 
jets in pA collisions at RHIC the beam quarks q* have small transverse momentum 
and the production of high-p quarks is not affected by virtual radiative corrections, see 
the discussion in Sect. VII B. Without further details we quote again first the dipole 
representation: 

2(2n) 2 da N (q* -> q(p )g) 



J d 2 rd 2 r'^(z q , r)^*(z q , r') exp[ip q (r - r') 



dz q d 2 p q 

{a 3 (r, r - z q r') + a 3 (r', r' - z q r) - a qq (r - r') - a qq (z q (r - r'))} , (HO) 

while the linear /c^-factorization result reads 
2(2n) 2 da N (q* - q(p q )g) 



dz q d 2 p q 



= J d 2 Kf(x A ,K) 



C 

{2^[Mz q ,p q )-^(z q ,p q + K)\ 2 

+ \*( Zq ,p q ) -V(z q , Pq + Z q K)\ 2 -\V(z q ,p q + K) -V(z q , Pq + Zq K)\ 
+ \V(z q ,p q + K) -*(z q ,p q + Z q K)\ 2 }. (Ill) 

A closer comparison of the p spectra of gluons and quarks from g* — > gq reveals interest- 
ing differences which hold in both the nonabelian and abelian, Ca — > 0, cases. Evidently, 
in the p-integrated spectra we would obtain the longitudinal spectrum of quarks simply 
by swapping z g — > 1 — z g = z q in the formula for the gluon spectrum and vice versa. The 
similar procedure for the p-spectra, as advocated in [21] would be entirely wrong. 



D. Quark spectrum from the breakup of the qg— Fock state of a fast quark on 
a nuclear target 

Routine application of the above described technique yields 

(2n) 2 da A (q* -> q(p q )g) _ 
dz q d 2 p q d 2 b 

AA 



S abs {b) J d 2 K<f){b,x A ,K){\q{z q ,p q ) -V{z q ,p q + K)\ 2 + \*(z q ,p q ) -*(z q ,p q + z q K)\ 2 } 

+ J d 2 K 1 d 2 K 2 <fr(b,X A ,K 1 )(fr(b,XA,K2)\^(z q ,p q + Ki) -^(z q ,p q + Z q K 2 )\ 2 . (112) 

As usual, the absorbed linear /cj_-factorization term is the counterpart of the large-A^ 
free-nucleon spectrum (111). Note also that the p ? -spectrum of slow quark jets from 
q* — > qg coincides with the spectrum (34) form g* — + QQ and the Cronin effect for the 
two cases will be identical: 

A qq (z q < V,x A ,P q ) = A Qg (z Q < l;x A ,p q ) . (113) 

For the numerical estimates for pAu collisions one must take Ac r onin(x A ,p q ) for j — 4 
in Figs. 8, 9. 

The spectrum of leading-quark jets, z q — > 1, equals 
(2n) 2 da A (q* -> q(p q )g) 



^=2S abs (b) J d 2 K(f)(b,x A ,K)\y(z q ,p q ) -V(z q ,p q + K) 



dz q d 2 p q d 2 b 

+ J d 2 K 1 d 2 K 2 <l)(b,X A ,H,i)(f)(b,X A ,K2)\^(z q ,p q + Ki) - ^(Z q ,p q + K 2 )\ 2 , 

(114) 

where the second term is similar to da^ of eq. (96) which emerged in the discussion of 
the spectrum of leading gluons from g* — > gg. Its antishadowing properties have been 
studied in section VII B and entail 

A qq (z q -> l;x A ,p q ) =A qq (z q < l;x A ,p q ) (115) 

The antishadowing effect for leading quarks from q ^ qg is one half of that for leading 
gluons from g ^ gg. 

The above derived ^-distribution quantifies the energy loss for the radiation of glu- 
ons, i.e., gives a solution of the LPM problem for finite transverse momentum p of the 
observed quark jet. The case of quark jets is important for its relevance to the large 
pseudorapidity, high-p hadron production in deuteron-gold collisions at RHIC. The nu- 
merical studies of the LPM effect contribution to nuclear quenching of forward high-p 
hadrons observed experimentally by the BRAHMS collaboration [18] will be reported 
elsewhere. 

45 



IX. DISCUSSION OF THE RESULTS 



Extending the technique developed earlier [2] we elucidated the pattern of k±- 
factorization breaking for the single gluon and quark jet production from g* — > QQ, g* — > 
99, 1* ~^ Q9 breakup on nuclear targets. In all the cases we derived an explicit nonlin- 
ear, from quadratic to cubic, /cx-factorization representations for the inclusive single-jet 
transverse and longitudinal momentum spectra in terms of the collective nuclear uninte- 
grated gluon density. The results are simple and elegant and have not been presented in 
this form in the existing literature on the subject. They also give a solution of the LPM 
problem for finite transverse momentum p of the observed secondary quark and gluon 
jet. One of the evident future applications will be to the nuclear quenching of forward 
high-p hadrons observed experimentally by the BRAHMS collaboration [18]. 

The unintegrated collective nuclear gluon density which defines the nuclear observ- 
ables is shown to change from the quark to gluon observables even within the same 
process what exemplifies the point that an ultrarelativistic nucleus can not be described 
by a single gluon density, the collective nuclear density must rather be described by a 
density matrix in color space [2]. The often used linear /cj_-factorization description of 
particle production off nuclei (e.g. [23], for the recent applications and more references 
see also [24, 25, 29, 30]) is definitely not borne out by the multiple scattering theory. 
There is a striking contrast to the single-quark jet transverse momentum spectrum from 
7* — > qq breakup on a nucleus for which linear /c_i_-factorization is fulfilled, and our anal- 
ysis shows that that the linear /cj_-factorization in 7* — > qq is rather an exception due 
to the color-singlet nature of the photon by which the single-jet spectrum accidentally 
becomes an abelian problem [2]. 

Although the linear A;_i_-factorization for hard processes in nuclei does not exist, for all 
breakup processes we identified the 'abelian' regime in which the transverse momentum 
spectra of soft quark and gluons jets with z q>g <C 1 take the process-dependent linear k±- 
factorization form. Specifically, the production of quark jets is described by the collective 
gluon density (f>(b,XA, «0 defined for the triplet-antitriplet color dipole. The production 
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of gluon jets is described in terms of the collective gluon density (fi gg (b, xa,k) defined 
for the octet-octet color dipole. The color-dipole representation for such a slow gluon 
limit coincides with that derived before for the slow gluon production by Kovchegov 
and Mueller [20]. The reason behind the recovery of the Kovchegov-Mueller form for 
slow jets is that in this limit the single-parton spectrum becomes an abelian problem 
[2]. Such a form is of limited practical significance, though. Evidently, the limit of 
z q <C 1 is irrelevant to the open charm production, in the gluon jet production too the 
dominant contribution to the observed jet cross section will come from finite z to be 
described by our nonlinear /c^-factorization. The full z-dependence is also mandatory 
for the evaluation of the LPM quenching of forward high-p jets. The production of 
gluons from g — > g ± g 2 has been discussed also in [22], their equation for the nuclear 
spectrum contradicts our general result (17). The earlier work on q — > qg by [21] starts 
with the two-parton inclusive spectrum. When compared to our general equation (17), 
in ref. [21] the coupled-channel S-matrix Sgf gg has illegitimately been substituted by 
a single-channel Glauber exponential. Despite this inconsistency, the final color dipole 
representation of Ref. [21] for the single-parton spectrum agrees with ours. 

Recently, the production of heavy quarks from the excitation q — > q'g, g — > QQ has 
been discussed in the color dipole framework in ref. [31]. We have however not been 
able to reconcile the results given there with our findings on the cancellation of spectator 
interactions. Very recent works [29, 32] start with the reference to the color glass con- 
densate model, but in actual applications the Glauber-Gribov multiple scattering theory 
is used and a quantity similar to the collective unintegrated glue of [2, 12] is introduced. 
Although no explicit formulas corresponding to eq.(31) for the heavy quark cross section 
are given, the authors are aware in part of our criticism [2] of the applicability of linear 
k± factorization in the saturation regime. 

Regarding the applicability domain of our formalism the considered processes corre- 
spond to production of parton pairs with Xb, c %a ~ 10~ 2 , where Xb, c are defined with 
respect to the nucleus. In the pA collisions at RHIC this requires that both partons were 
produced in the proton hemisphere. The consistent discussion of higher order processes 
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and of the emerging BFKL evolution of nuclear spectra remains an open problem. Gen- 
erally speaking one will have to address a situation in which the coherence of more than 
two partons over the nucleus is important, and the question is in how far the rescatter- 
ing of partons faster/slower than the measured jet is reabsorbed into the (non-)linear 
evolution of the beam/target effective gluon densities. Such evolution effects would per- 
haps not feature too prominently at RHIC energies, where the rapidity span overlaps 
well with the applicability region of our approach. Still such studies will be important 
to establish the viability of our approach to higher energies such as LHC. They would 
also sharpen our understanding of the similarities and differences of our color dipole for- 
malism and the Color Glass Condensate approach (for recent explorations of the energy 
dependence of the Cronin effect for slow gluon production within those models, see e.g. 
[24-27, 32-34]). 

The most obvious application of the derived nonlinear nuclear /c_|_-factorization is to 
nuclear effects in the single-jet spectrum which are usually referred to as the Cronin 
effect. Our derivation unambiguously relates it to the antishadowing property of the 
collective nuclear gluon density [2, 12]. The Cronin effect was shown to exhibit interesting 
variations across the phase space. Based on the numerical studies of the antishadowing 
effect in the collective nuclear gluon density [12] we presented numerical estimates for 
the Cronin effect. One feature of these estimates is noteworthy: The realistic models 
for the unintegrated gluon density of the proton [14] suggest a somewhat small value of 
the saturation scale, Qa ~ 1 GeV, even for the heaviest nuclei. Still, the Cronin effect 
is found to reach its peak value at |p| ~ (2 4- 2.5) GeV, which is well within the pQCD 
domain. 

To summarize, our analysis establishes nonlinear /cj_-factorization as a universal fea- 
ture of hard production off nuclear targets. The presented formalism can readily be 
extended to the nuclear dependence of jet-jet correlations and nuclear quenching of for- 
ward quark jets, the corresponding work is in progress. 

Acknowledgements: We are grateful to B.G. Zakharov for illuminating discussions 
and suggestions. 
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Appendix A: The derivation of 



Here we expose some of the technicalities behind the derivation of the master equa- 
tion (10) and the reduction of the four-parton ^-matrix to the two-body ^-matrices in 
equation (14). The color-dipole properties of excitation amplitudes are elucidated and 
the calculation of the S'-matrices is greatly simplified in the qq representation of the 
gluon interactions. We demonstrate the principal virtues of this technique on an example 
of g* — > QQ on a free-nucleon target, for the extension to a nuclear target see [2]. In 
the color space the g a QQ vertex contains the SU(N C ) generator t a , so that in the quark 
color space the contribution of the diagram of fig. 3b is proportional to S (b b )t a (b c ) , 
whereas the contribution of the diagram of fig. 3c is proportional to —t d Sd a (b). Here 
S^(b c ) describes the interaction of the antiquark, and S da (b) = (gd\S g (b)\g a ) describes 
the transition between the two gluons in color states a and d. Consequently, the properly 
normalized excitation operator is 

M a (b b , b c , b) = -j±==(s(b b )t a S\b c ) - t d S da (b)) (A.l) 

where CfN c = Trt a t a . The above vertex operators t a have been suppressed in the 
expansion (6). 

Now we note that 

S da (b) = ^Tr (t d S(b)t a S\b)) (A.2) 
and, upon the Fierz transformation, 

t d S da (b) = S(b)t a S^(b), (A.3) 

which leads to a particularly simple form of the excitation operator: 

M a (b b ,b c ,b) = -^==(s(b b )t a S\b c ) - S(b)t a S\b)) (A.4) 

Note how M a (b b , b c , b) vanishes if b b = b c = b: in this limit the QQ behaves as a 
pointlike color-octet state indistinguishable from the gluon and its QQ structure can not 
be resolved. What enters the integrand of (10) is Tr (^Ml(b' h , b' c , b')M a (b b , b c , bfj. 
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Before proceeding further we notice that the S'-matrix for the color singlet qq dipole 
equals 

S 2 (h - b c ) = ^Tr (s{b b )S\b c )) = 1 - F 2 (B, b b - b c ) , (A.5) 

where r = b b — b c is the relevant dipole size and T 2 (B,b b — b c ) is the corresponding 
profile function. Upon the integration over the dipole impact parameter, B, one finds 
the color-dipole cross section 

a g g(x, r) = 2 J d 2 BT 2 (B, r) . (A.6) 

The excitation on a free nucleon target is described to the single-gluon exchange approx- 
imation. In this approximation 

Sl(b b -b c ) = l- 2T 2 (B, b b - b c ) , (A.7) 

Now we can identify the in the integrand of (10) and give explicit expressions in 
terms of T 2 (B, r). Making use of 

Tr (t a At a B) = -(Tr (a) • Tr (b) - -J-1Y (AB)) (A.8) 
2 N c 

and applying the unitarity condition (13) wherever appropriate one readily finds 

S< 2 \b,V) = -^Tr(S(b')t a S\b')S(b)t a S\b)) = l-^T 2 (B,b'-b) (A.9) 

S (3 \b b , b c , b') = ^Tr (S{b b )t a S\b c )S{b')t a S\b')) (A.10) 

= 1 " |^(r 2 (B,6'-6 6 ) + r 2 (B,6'-6 c ))-^^r 2 (B,6 b -6 c )| 

S {3 \K,b' c ,b) = ^-L-Tr (s(b> b )t a S\b> c )S(b)t°S\b)) 

= i - [^{UB,b-b' b ) + r 2 (B,b-b'j)-^^r 2 (BX-K)Y 

(A.ll) 

In the one readily identifies the profile functions of the qqg cross section as introduced 
in [37]. The case of is more tricky: 

S^(b b ,b c ,b' b ,b' c ) = ^-Tr (S(b b )t a tf(b c )S(b' c )t a S\b> b )) 
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= ^rriKt 1 - b c - K) ~ ^(b, b b - b' b )) 

- ^Tr [S(b b )S\b c )S(b' c )S\b' b )]}. (A.12) 

The last term in (A.12) is an operator in the qq color space and in general is 
the coupled-channel operator derived in [2]. In the case of the single-particle particle 
spectrum b' c = b c and the straightforward application of the unitarity condition (13) 
yields 

S^(b b , b c , b' b , b c ) = S 2 (b b - b' b ) = 1 - r 2 (B, b b - b' b ) . (A. 13) 

As quoted above the results are for the free nucleon target, for the nuclear target one 
must apply the Glauber-Gribov exponentiation to the relevant S^ n \ 



Appendix B: 

In all the transverse momentum spectra the recurrent quantity is \ty(z,p) — *&(z,p — 
k)\ 2 . This quantity for the incident photon is found in [2, 5]. For transverse photons 
and flavor /, i.e., for the 7* — > QQ excitation, 



\V(z,p) -^(z, P -k)\ 2 = 2N c e 2 Q a en 

\ 2 



x<[z 2 + (l-*) 21 ' P P K 



2 



p2 _|_ £ 2 (jy _ ^2 _|_ £ 

+m 2 Q (^-^-- \= 2I 1 ' ( R1 ) 

Q \p 2 + e 2 (p-k) 2 + e 2 J J 

The term oc m 2 Q must only be kept for heavy quarks. The same quantity for the g* — > QQ 
is obtained from (B.l) by the substitution N c e 2 QOt em — > TpasiQ 2 ) an( i ^ ne substitution 
in (B.2)of the virtuality of the photon Q 2 by the virtuality of the beam gluon Q 2 ,. For 
the general case a — * bc 

e 2 = z b z c Ql + z b m 2 c + z c m\ , (B.2) 

where Q 2 is the virtuality of parton a and the mass dependent term is important for 
heavy flavours. Now note that the factor T F [z 2 + (1 — z) 2 ] which emerges in the first 
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term in the rhs of (B.l) is precisely the familiar splitting function Pq 9 (z). For all other 
cases \fy(z,p) — ^(z,p — k)\ 2 is obtained from that for g* — » QQ by the substitution 
of Pq 9 (z) by the real-emission part of the relevant splitting function for z < 1 found 
in all textbooks [1]. Take for instance the fragmentation of light quarks q — > qg. The 
\iff(z q ,p) — ty(z g ,p — n)\ 2 which enters the spectrum of quarks must be evaluated with 
the splitting function 

P qq (z q ) = C F ^- y (B.3) 

The \ty(zg,p) — ty(z g ,p — k)\ 2 which enters the spectrum of gluons in the same process 
must be evaluated with 

PM=C F 1 + {1 - Za)2 . (B.4) 

In the fragmentation of gluons g — > gg one must take \^(z,p) — *&(z,p — n)\ 2 with the 
splitting function 



1 — z z ' 

+ - + z(l-z) 

z 1 — z 



(B.5) 



If e 2 is negligible small compared to p 2 , then one can use the large-p approximation: 

(p p-K V K 2 

\P 2 {P — K) 2 J P 2 (P — ^) 2 ' 

and it is worth to recall the emerging exact factorization of longitudinal and transverse 
momentum dependencies which is a well known feature of the high energy limit. 
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